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, Abstract 

■ An exact solution of the energy shift in each quantum mechanical energy levels in a one dimensional 
symmetrical linear harmonic oscillator has been investigated. The solution we have used here is firstly 
derived by manipulating Schrodinger differential equation to be confluent hypergeometric differential 
equation. The final exact numerical results of the energy shifts are then found by calculating the final 
analytical solution of the confluent hypergeometric equation with the use of a software (Mathcad Plus 
6.0) or a program programmed by using Turbo Pascal 7.0. We find that the results of the energy shift 

' in our exact solution method are almost the same as that in Barton et. al. approximation method. 

Thus, the approximation constants appeared in Barton et. al. method can also be calculated by 

■ using the results of the exact method. 

> 

On 

' 1. Introduction 

o ■ 

In 1990, Barton et. al. had already published an approximation method to calculate the energy shift in each quantum 
&\ mechanical energy levels in a one dimensional symmetrical linear harmonic oscillator. They found a formula of the 
ON ' energy shifts either for the ground state or for the excited states calculated in a finite interval L. The explanation 
about the use of their approximation formula is for bound states in one dimensional symmetrical linear harmonic 
O , ; oscillator with a potential V (x) which is proportional to x 2 . W However promising their method, the exact energy 
shifts in each energy levels have not been solved until now yet. In this paper We use another method to solve 
the problem. Instead of using the Barton et. al. approximation method, We apply a new method by changing 
Schrodinger differential equation to a confluent hypergeometric differential equation. After finding the confluent 
^ ( - hypergeometric differential equation, We solve exactly the equation and then use a certain program or software to 
calculate the final exact energy shifts. 

This paper is organized as follows. In section 2, we will perform the exact solution method of the energy shift 
in each quantum mechanical energy levels in a one dimensional symmetrical linear harmonic oscillator and will also 
compare our results with the results in ref.[l]. Section 3 is devoted for discussions and conclusions. 



2. Exact solution method of the energy shift in each quantum mechanical energy levels in a one dimensional 
symmetrical linear harmonic oscillator 

A particle with mass m moving in a one dimensional space from — oo to oo with the influence of a potential V (x) obeys 
the following Schrodinger differential equation! 1 ! 

+V(x)\ipo(x) =E i/)o(x), -oo<x<oo. (2.1) 

On the other hand, if the motion of the particle is restricted by a finite interval — ^ < x < h in the same potential, 



*email: hendry202@cyberlib.itb.ac.id 



1 



the wave function tp (x) and the energy E can be obtained by solving the following time independent Schrodingcr 
differential equation 



2m dx 2 



+ V(x) )ip(x) = Eip(x). 



L L 

2 <X< 2 



(2.2) 



Even though the energy shift AE = E— E is small for the large L, this case is not included in a classical perturbation 
theory. 

Based on this simple physical phenomena, Barton et. al. derived generally an approximation method for calculating 
the energy changes or the energy shifts in each quantum mechanical energy levels including the systems in linear 
harmonic oscillator and hydrogen atom. They got the following formulas 
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, for ground state, 



and 
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for excited states. However, for the calculation of the energy shift in each quantum mechanical energy levels in a one 
dimensional symmetrical linear harmonic oscillator with a potential V (x) = muj 2 x 2 /2 } they found 
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, for ground state 



(2.4a) 
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where V'o ( x ) — ( ^1/^1/2 ) ex P ( — §p ) , I — (^j) 1 an d the initial energy is E =^?kv. And for the excited states, 
they got 



E, 



(n) 



(2n+ l)7r 1 /2 n !2« ^ I 
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(2.4b) 



where (x) = ( 7 r 1 / 2 /2 2 n!) 1/2 H n (f ) exp and £< n) = (§ + n) huj, n = 1, 2, 3, ... . For the large x, H n (f) ~ 

2 n (f)" 

Now, we are going to change the Schrodingcr differential equation to become a confluent hypergeometric differential 
equation. We provide eigen function Hip = Eip for a one dimensional symmetrical linear harmonic oscillator with a 
potential V (x) = muj 2 x 2 /2 into two areas -\ 2 ^ 



and 



2m dx 2 



h 2 d 2 



mu! + x I ip+ (x) = Eip+ (x) , forx > 0, 



'2~mA~x 2 + 2 mLU - x2 ) ^~ ^ = ^ ' imX < °' 



By substituting parameters = A, and ^ =aora= p, we find 



dx 2 



- [A - a 2 x 2 } iP± = 0, for ± x > 0. 



(2.5a) 



(2.5b) 



(2.6) 
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We then introduce parameter A = a (2n' + 1) and a free variable z = Xyfa to eq.(2.6), we get 



d 2 ip n ' (f) 
dz 2 



[2ri + 1 - z 2 ] ij) n , (z) = 0, for ± z > 0, 



and 



E n . =hcj[-+n / ) . 



It means that 

o->_ oj 2n' + 1 



u> 4 



2n' + l 



r = 1, 



(2.7) 



(2.8) 



(2.9) 



which states that our solution is for a one dimensional symmetrical linear harmonic oscillator. 

If ip n > (2) is the solution of eq.(2.7), then ip n i (—z) is also the solution. On the other hand, if the Wronskian 
W (ip n i (xy/a) , ipn' (~ X V&)) is not the same as zero, then both solutions can be used as a basic solution for imple- 
menting boundary conditions whether in x — > ±-| or in x = 0. To solve eq.(2.7), we firstly substitute the following 
function 

Vv (z) = exp 



H n , (z) , 



we then get 



dz 2 



-2z( 



\dz 



2n' 



H n . (z) = 0, 



(2.10) 



(2.11) 



where H n > (z) is a Hermite function. By introducing a new variable t — z 2 , we find a confluent hypergeometric (CH) 
differential equation 



dz 2 



1 



- t 



H n , (t) = 0, 



which has a general solution in a linear combination of the following two CH functions! 3 ! 



B 



(w) 



1F1 



1-n' 3 



In eq.(2.10), the Hermite function i?„/ (z) is just a special function of the following equation 



«,(.)- 2 "' r(l; 
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Fll ~T ; 2 ;z J + r(-£) 



r(|(i-n')) 

Here, we have chosen constants A and B as follows ^ 

2«'r(i) 



l-n' 3 



2 5 2 ;Z 



A = 



and 



B = 



r(i(i -«'))' 



2 n 'r(-|) 



r(-£) ^ 

The boundary conditions in 2; are 



(2.12) 



(2.13) 



(2.14) 



(2.15a) 



(2.15b) 
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H n , (0) 



2"'r (|) 
r(i(i-n')) 



and 



(Q) = 2"'r(-i) = 
dz r(-f) r(-f) • 

We can also show that! 3 ! 



W (ip n > (x\fa) , ip„> (-xy/aj) 



2 n 



r(-n') 



- exp (z 2 ) . 



(2.16a) 



(2.16b) 



(2.17) 



Hence, if n' is not exactly the same as n' = 1, 2, 3, ... , then both of the solutions H n i (±z) have a linear independent 
property. If this is work, then both of the solutions can be chosen as a basic solution of eq.(2.7), we get the general 
solution as 



(z) = [AH n . (z) + BH n , (-*)] exp 



1 



Here, the boundary condition for the eigen function in ±z = ±^yfa is 
L 



i> n , (±z) = il> n , \± -V^J = 0. 
By substituting the boundary condition in eq.(2.19) into eq.(2.18), we find 



AH n , 



'a 



BH n > 



exp 



a = 0. 



Eq.(2.20) can be divided into two equations with two conditions as follows 
(1). For the condition with even parity: 



Vv 



L 



a = 



AH n , 



L 



BH„i 



L 



exp 



0, 



or 



2"'r(|) 

r(i(l-n')) 
(2). For the condition with odd parity: 
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= 0. 



In both eq.(2.21a) and (2.21b), hermite functions H n > (±-| v / a)are 
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and 



H r A--^a) = 



2"'r(l) / „' i L 2 
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where 
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From eq.(2.21b), we find a solution for the odd parity, 

n' = 1 + 2m'. 

On the other hand, the even parity is found by calculating eq. (2.21a), we then get 
n' = 2m'. 



(2.23b) 



(2.24a) 



(2.24b) 



Here (in eq. (2.24a) and (2.24b)), n'and m'are positive real values due to the boundary condition z — ±^y/a where 
ip n i (±z) = vanishes in z = ±^\fa. n'in eq. (2.21a) and (2.21b) can exactly be calculated by using mathcad plus 
6.0 or a program programmed by using Turbo Pascal 7.0. The results of the calculation are provided in table 1 and 
2. 

The requirement in eq. (2.24b) produces the bound states quantum mechanical energy with even parity, 



E n ' 



Two, n = 2m . 



(2.25a) 



and the requirement in eq. (2.24a) produces the bound states quantum mechanical energy with odd parity as follows 
I s 



fko, nl 



2m' . 



(2.25b) 



On the other hand, for the calculation of the quantum mechanical energy in the boundary z — ±oo, and the 
boundary condition i[> n (±z) = finite, we can use an asymtotic formula of H n (±z) according to ref.[3] as follows 



37T 

H n (±z) ~ (±2z) n , for z — > ±oo or \z\ — > oo if |arg(z)| < — , 



(2.26a) 



and 



H n (±z) ~ (V^)exp(z 2 )(|z|)- 



for 



±00. 



(2.26b) 
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By substituting the conditions in eq.(2.26a) and (2.26b) into eq.(2.18), we find the bound states quantum mechanical 
energy with even parity, 

E n = (n+^jhu, n = 2m = 0,2,4,6... , (2.27a) 

and for the odd parity, we get 

E n =(n+^jfiu, n = 1 + 2m = 1,3,5, ... . (2.27b) 

Hence, for a one dimensional symmetrical linear harmonic oscillator in the boundary z = ±oo , we get its wave 
function as 

ip n (x) = [A n H n (x^/a)] cxp (-\ax 2 ^j , (2.28) 

where the constant parameter A n can be found by a normalization. 

Now, the energy shift in each quantum mechanical energy levels in the one dimensional symmetrical linear harmonic 
oscillator can exactly be calculated by formulating a simple formula related to eq.(2.25) and eq.(2.27) as follows 

A£ (n) = E n , - E n 

= (ri - n) Tvuj, (2.29) 
where n — 0, 1, 2, 3, 4, 5, ...(positive integer) and the values of parameter n' depend on jy/a, for example : for 

n' = 1.55xl0~ 15 , (l + 1.082xl0~ 13 ) , (2 + 3.671xl0" 12 ) , (3 + 0.805xl(T 10 ) ,etc. 
Some results of the energy shift are attached in tabic 1 and 2. 



k = \y/0L 


Barton ct. al. results (1990) 
(see : eq.(2.4a), AE^) 


Our results 
(A£(°) = E , - E ) 


1 


0.415(1 + (i)) Too 


0.798 Two 


3 


4.177xl0- 4 (l + 0(±)) Two 


3.911xl0- 4 Two 


4 


5.079xl0- 7 (l + O( ? i f )) Two 


4.908xl0~ 7 Two 


5 


7.835xl0- ii (l + O( T 4^)) Two 


7.671xl0- li ?iw 


6 


1.570xl0- ib (l + O( T i r )) Two 


1.550xl0~ 15 Two 


7 


4.141xl0- 2i (l + O( T |g)) Tiuo 


4.098xl0- 2i ^u; 


10 


4.197x10-^(1 + 0(^t)) Two 


36.769xl0- 43 ?iw 



Tabic 1. The energy shift in the ground state of the one dimensional symmetrical linear harmonic oscillator. 
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k = ±y/a 


Barton ct. al. results (1990) 
(see : eq.(2.46), AE^) 


Our results 

(A£7(») = E n , - E n ) 


3 


AE^ = 7. 52xl0~ 3 A huo 
AE^ = 6.767xl0~ 2 ATiuj 
where A = (l + 0(^)) 


AE^ = 6.081xl0- 3 Two 
AE^ = 4.119xl0- 2 ?iw 
we get A — 0.606 


5 


AE^ =3.918xlO- 9 Bftw 
AE^ =9.794xlO- s Bhuj 
AE^ = 1.632xl0~ 6 Bhco 
where B = (l + 0{± 5 )) 


A£W =3.672xl0- 9 ?iw 
=8.402xl0- 8 ?iw 
AEW =1.221xl0- 6 7iw 
we get B = 0.848 


6 


AE {1 ) =1.13x10- 13 C^j 
AEW =4.07x10- 12 C^ 
A£( 3 ) =0.98xlO- 10 C^ 
where C = (l + (^)) 


AEW =1.08xl0- 13 fiw 
AEW =3.67xl0- 12 ?ia; 
AEW =0.80xl0- 10 ?iw 
we get C =0.892 



Table 2. The energy shift in the excited states of the one dimensional symmetrical linear harmonic oscillator. 



3. Discussions and Conclusions 

We have presented an exact solution method of the energy shift in each quantum mechanical energy levels in a one 
dimensional symmetrical linear harmonic oscillator by using the solution of the confluent hypergeometric differential 
equation. We can conclude that the relationship of the energy shifts is as follows : 

AE {0) < AE (1) < AE {2) < AE (3) < ... < AE {n - 1] < AE (n) . 

According to the comparison between our results and Barton ct. al. results shown in table 1 and 2, we get that 
their results are almost the same as our results. However, our results are the exact results. On the other hand, 
our results will be more accurate than their method if in our calculation, we add the series of gamma function in 
eq. (2.23a) and (2. 23b). The longer gamma function series, the more accurate our results. Based on the comparison, 
we also get the approximation constants (obliquity factors) appeared in Barton et. al. method. 

In conclusion, our method is work without a detail calculation of the wave function in each energy levels. 
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